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Page 7: The Important derivatives and are obtained 

incorrectly. The corrected values for these derivatives should read 
as follows : 


xk E^fBO) - M^BC) 
2 M2 _ 1 


(4a) 


, _ jtA B^CeC) -M 2 h(BC) 

16 M2 - 1 


(5a,) 


The same correction should be made for these derivatives in table I 
in the column "Principal, body axes." In the column. “Stability axes" 
of table I the corrected expressions are 



* A E» , (BC) - M2h(BC) 
2 M2 _ 1 


u. = asifi + s 5gV'.(S°>. 

^ 16 \ c / M 2 — 1 


The error in the derivation consists in the assumption, carried 
over from reference 1, that the surface potential, equation ( 3 ) j is 
not altered by a small normal acceleration. This assumption is true 
for the narrow triangles treated in reference 1 but fails for the 
general triangles treated, in the present paper. The assumption is 


equivalent to the neglect of the time dependency' terms 


2 7 
a 2 




and ~ 0^ (where t represents time and a is the speed of sound) 
& 

in the linearized partial differential equation for unsteady motions: 
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- 0 - 


yy 


t 


zz 


+ 5 0xt + i 0tt = 0 


(a) 


which, leaves only the steady-state (Prandtl-G-lauert) equation 


B 2 ^ -d>-$ =0 (b) 

~yy ”zz v 

The correct potential to replace equation ( 3 ) must satisfy both 
equation (a) and the boundary condition on the wing 



(c) 


In an unpublished paper. Me. Clifford S. Gardner has, in effect, shown 
that a suitable solution is 


where If is the steady— state potential corresponding to a unit 
pitching velocity about the y-axis and X i B the Bteady-etate 
potential corresponding to unit angle of attack. Both \|r and X 
satisfy equations of the form (b). That equation (d) is a solution 
can be verified by direct substitution into equations (a) and (c). 
Thus, Gardner has shown that the time— dependent potential for an angle 
of attack &t may be compounded of two time— free, or steady-state, 
potentials, one for a constant- angle of attack and the other for 
steady pitching. 


The lift distribution at time t = 0 for the angle of- attack at 
is obtained from the potential by 

AT = 2 P (V0 X + 0 t ) 



(e) 
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where 

(AP) lift distribution, for unit pitching velocity 

(AP)^ lift distribution for unit engla of attack 

AP* Talue of AP per unit at, used in deriving the 

incorrect equations ( 4 ) end (5) 

Integration of equation (c) over the plan form to obtain the 1 lf t 
moment and reduction to coefficient form yields 


**» - £ % + 8 £ 'V. - 


(f) 


r . _ Mf 2 M. 

- B 2 % + e %s2 


u 


n 


dx dy - -4- 


plan 

form 




B2 


<ss 


(g) 


0=1 


where the * designates the incorrect values in equations (4) end (5) 
respectively 3 and the y-axis is taken through the center of gravity. 
Values of Ct , CL, s and C m , may be obtained from table I 
\* ^ ' x * y. 

and. ^C j j and ^ , from equations ( 4 ) and (5) . The 


quantity 


( - — -) is obtained by setting a = 1 and a = c($ c + x) 

[ipyy V3 ' 

X2 ' cx=l 

in equation (l) . Substitution and integration then yields the 
corrected values for and (equations (4a) and (5a)) as set 

forth at the beginning of this erratum. 
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Page lh: The section entitled "DERIVATIVES Cy^, Cng, CY r , ari ^ < ~' r \ L . ” 

and the part of table I in which these derivatives are presented have 
been found to be in error. The null result deduced for these 
derivatives (in the neglect of friction drag) is incorrect as a result 
of the omission of a compressibility factor in the treatment of the 
leading-edge suction force. Thus, on page l4 (last sentence of 
first paragraph), the remark that the initially symmetric distribution 
of leading-edge suction persists in sideslip (for Bmall values of 
sideslip) is a misstatement. The initially .symmetric distribution of 
leading-edge velocity persists in sideslip, but compressibility upsets 
the symmetry of the leading-edge suction . According to equation (17) 
the suction per unit length of edge is 


f = *pK 2 \Jl - M’ 2 


For infinitesimal sideslip the constant E, related to the edge 
velocity, is unchanged, but M’ , the component Mach number perpendicular 
to the edge, is altered: M* increases on the right edge and decreases 

on the left edge. The err oh -under discussion consists in the neglect 
of the change in M* . 

Because of the change of M* with sideslip angle 0 the edge 
suction may be written 


f 



( 21 ) 
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By differentiation of equation (17), with K constant and 
M* = M sin(<= ± 0), there results 


s-et 


where the upper part of the double sign refers to the right- edge and 
the lower part to the left edge. 

The quantity fp_o is obtained by setting p = 0 in the last 
equation of page 12: 


■P-o 


TtpCxc^V 2 \Jl - B^C 2 

2 §C f (BC}] 2 


(23) 


Substitution of equation (23) in the last term of equation (22) and 
simplification, with tan e = C, yields 


f = f 


0=0 


T P 


7 rpY 2 a g 3:C g M g 


2[e*(BC)J 2 \l 1 - B^C 2 


(2 4) 


Equation (24) gives the suction per unit length of edge for a 
triangular' wing with an angle of sideslip 0. 

For the case of a small angular -velocity of yawing r the edge 
suction may be approximated by 


f ’ + r (£L 


where f r= o I® same as fp=o — and is given by equation ( 23 ) . If the 
center of rotation is (^°> 0, 0^ , the component Mach number normal to the 
edges is — * r - ■»•••*■ 


M* = M^sin 6 ± 7 ( 3 0 cos 6—2 sec 


it y 
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This value of M* is to he incorporated in equation ( 17 ) for f 

before carrying out the indicated differentiation — . The final result is 

5r 

fi-c — 1 sec 2 e') irpV a 2 xC 2 M 2 

f = f j3=0 Tr - 4 ■— == (25) 

2[e’ (BO)"] 2 \|1 — B 2 C 2 

The suction equations (24) and ( 25 ) have been used in procedures 
similar to those presented on page 13 to yield corrected values of 
the stability derivatives CL. , C , CL. , and a correction term to be 

*13 Hp x r 

added to the right-hand side of equation (20). The results are 
incorporated in a correction to table I in the column labeled 
"Principal body axes". Since one or another of these four derivatives 
figures in their transformation to stability axes, several of the 
equations in the column of table I labeled "Stability axes" ha ve also 
required correction. 

The yawing derivatives Oy and as corrected herein are 

st.-m subject to the uncertainty inherent in the limitations of the 
linearized theory for yawing motion pointed out on pages 9 and 10. 


Page 18: Formulas 

changed to read as 


for G n p * C np> °Yp* ^ % 

shown in attached correction to table I 


should be 



Corrected version of part of table I following the derivative C 7 . 


Source 
(if not 
derived 
herein) 



Principal body axes 
, ^origin at 


■rfr a2A2^Q(BC) 

r 


Stability Axes 


/ p \ 

(origin at distance x Q g ahead pointy 


J^Udc) + (je + 


- + ib) 


J(BC) 


~ Ms + ife 


§x + it + i ^) J(BC) - is I(BC) ] - + ^ 


- + * re + & ¥) &,,(bc) - 


- + §^' BC > - 


2 A Kbc) - + A + 2A! +A ^c S+ 8^ 


nar ^ I(BC) - ^ ^ 

32 9 \A 8 256 


c Ac 


Q(BC) 


- J a 2 AJ^Q(BC) 


^ J(BC) 


4- a 2 A^J®Q(BC) 

2V 


- f a 2 A^Q(BC) 
4 


J(BC) 


J(BC) +(~ + ^_ 
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TECHNICAL NOTE NO. 1572 


STABILITY DERIVATIVES OF TRIANGULAR WINGS AT 
SUPERSONIC SPEEDS 

By Herbert S. Ribner and. Frank S* Malvestuto, Jr. 


SUMMARY 


The analysis of the stability derivatives of low— aspect— ratio 
r5 angular wings at subsonic and supersonic speeds, given in NACA TN 
'o. 114-23, Is extended to apply to triangular wings having large vertex 
ogles and traveling at supersonic speeds. The lift, rolling moment 
je to sideslip, and damping in roll and pitch for this more general 
case have been treated elsewhere on the basis of the theory of small 
disturbances. The surface potentials for angle of attack and rolling 
taken therefrom are used to obtain the lift due to downward acceleration, 
the several side— force and yawing— moment derivatives that depend on 
leading-edge suction, and a tentative value for the rolling moment due 
to yawing. All the known stability derivatives of the triangular wing 
at sxiper sonic speeds, regardless of source, are summarized for convenience 
and presented with respect to both body axes and stability axes. The 
results are limited to Mach numbers for which the triangular wing is 
contained within the Mach cone from its vertex. The spanwise variation 
of Mach number in the case of yawing is neglected, although the effect 
must be of importance. 


INTRODUCTION 


An earlier investigation (reference 1) has provided theoretical 
stability derivatives of low-aspect-ratio wings of triangular plan form 
at subsonic and supersonic Bpeeds. The restriction to low aspect ratio 
was a consequence of the limitations of the theory. Several investigators 
have si »ce obtained pressure distributions for angle of attack, rolling, 
itchin,. , and sideslip at supersonic speeds (references 2 to 6 and 
lpublir , ' i ed analyses ) , without restriction to low aspect ratio. These 
■>rivati s have employed variants of the linear theory of supersonic 
ow and ave, in fact, constituted important steps in the development 
the i ory. 

If the rotations are taken about the vertex, the pressure distribu— 
on for each motion in the more general case is found to have the same 
,hape as the corresponding low— aspect— ratio approximation, so long as 
the triangular wing is contained within the Mach cone from the ve rtex. 
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The magnitudes differ by factors which are f tract ions solely of the ratio 
of the tangent of the semivertex angle of the triangle to the tangent of 
the Mach angle. The same similarity- exists between the distributions of 
surface potential. It is thus relatively simple to extend most of the 
derivations of reference 1 to remove the restriction to low aspect ratio 
for supersonic speeds. Such an extension is made in the present paper. 

The lift-curve slope, the damping in roll and pitch, and (in effect) 
the rolling moment due to sideslip have been evaluated in references 2 
to 6, so that the principal contributions of the present paper are the 
several, side— force and yawing-moment derivatives and a tentative value 
of the rolling moment due to yawing. All the known stability derivatives 
of the triangular wing at supersonic speeds, regardless of source, are 
collected herein for convenience and presented with respect to both body 
axes and stability axes. Wings with dihedral are not treated (although 
they were included in reference l) and the results are limited to Mach 
numbers for which the wing is contained within the Mach cone from its 
vertex. 


SYMBOLS 


x,y,z 

u,v,w 


V*<L,r 


Y 

M 

M* 


B 

a 

3 

€ 


M- 


rectangular coordinates (fig. l) 

incremental flight velocities along x— , y— , and z-axes, 
respectively, figure 2j Induced flow velocities along 
x—, y— , and z-axes of figure 1, respectively 

angular velocities about x— , y— , and z— axes, respectively, 
figure 2 


flight speed 

stream Mach number (V/Speed of sound) 
component Mach number normal to wing leading edge 

cotangent of Mach angle (n/m 2 — l) 
angle of attack (Flight w/Y) 
angle of sideslip (Flight v/v) 
semivertex angle of triangle 
Mach angle (cot”’ 1 ' \ / M 2 — l) 
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AP 

P 

a 

b 

c 

c 


local pressure difference between lower and upper surfaces of 
airfoil, positive In sense of a lift 

density of air 

semiwidth of triangle at distance x from vertex 
span (base of triangle) 
root chord (height of triangle) 


S 

0 

TJ = COS' 
* - 

E* (BC) 

F»(BC) 

e"(bc) = 

G(BC) = 


mean aerodynamic chord (is- / 

\ s U C 

edg* slops (| - g - $ - |f) 

aspect ratio Qrty 

area of triangle ^ibc) 
velocity potential 


b/2 


Q 

(Local chord) dy 


•) 


-1 Z 
a 


- B 2 C 2 


complete elliptic^ integral of the second kind with 
modulus k: 


; ^ 2 \Jl — k^sin 2 ! dz^ 


complete elliptic __ integral of the first kind with 

JC 

modulus kj ( / 2 


dz 


0 


E'(BC) 


1 - B^C 2 


(1 - 2B 2 C 2 )E*(BC) + B 2 C 2 F*(BC) 


H(BC) = 3G(BC) - 2E"(BC) 


CvjIoo 



k 
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i(bc) = 


2(1 - B 2 ^ 2 ) 


(a - b 2 c2)e»(bc) - bSc^f^bc) 

J(BC) = E"(BC)l(BC)v v /l -B 2 C S 

c 'Dust ant defined in equation (l6) . 
yawing moment 
lateral force 


m 


'n 


T N 

s 


suction force per unit length of edge 
, Lift ' 

lift coefficient 


r dV^S 


pitching-moment coefficient 


roll ing-moment coefficient 


yawing-moment coefficient 


lateral-force coefficient 


prof lie-drag coefficient 


induced surface velocity normal to wing leading edge 
perpendicular distance of point (x,y) from wing leading edge 



eg 

Subscripts 


distance of center of gravity forward of 


B 

L 


right edge 
left edge 
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Whenever <x, a, q., p, g, and r are used as subscripts, a 
nondimens lonal derivative is indicated and this derivative is the slope 
through, zero. For example. 



A dot above a symbol denotes differentiation with respect to time. 
All angles are measured in radians. 


ANALYSIS 

SCORE 


The stability derivatives of triangular wings at supersonic speed 
that have been treated theoretically herein or elsewhere are listed in 
table I, together with the expressions that have been found for them. 

All the derivations male toss of body axes. The derivations that follow 
give the values -with reference to the principal body axes of figure 2 
with origin at the aerodynam i c center f§c, Q»Oy, Conversion has been 

made to the system of stability axes shown in figure 3 with origin a 
distance x C g ahead of the point. Table I comprises parallel columns 


which present formulas relative to both systems. The expressions are 
limited to Mach numbers for which the triangle is contained within the 
Mach cone from its vertex. 


DERIVATIVES C T , C T . , and C„ 
% L a “a 


The pressure distribution on a thin delta wing at an angle of attack 
in a supersonic stream has been obtained in references 2 to ^ by the 
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linearized, theory without restriction on the vertex angle of the triangle. 
The approximation originally given for the slender (low-aspect— ratio) 
triangle (reference 7 ) and used as the "basis for reference 1 is found to 
apply to the general case upon division "by a constant (an elliptic integral) 
that depends on the ratio of the semivertex angle to the Mach angle. That 
is. 


AP UaCo, n) 

| pT 2 E* (BC ) - y2 

where E*(BC) is the complete elliptic integral of the second kind with 
modulus 

k = \Jl - B^ 2 
% 

= \Jl - (m 2 - ^ C 2 

Thus the lift-curve slope for the more general case is the value given by 
references 7 and 1 divided by E*(BC): 


c L 

a 2E , (BC) 

= | AE” (BC) (2) 


The surface potential given in equation ( 3 ) of reference 1 is likewise 
extended to include nonslender .triangles at supersonic speeds upon division 
by E*(BC'). The revised potential is 


„ * s u aj 

E'(BO) 


_ + 


Va 


x/a 2 — y 2 


E* (BC) 

pends o 

(ratio of the tangent of the semi vertex angle of the delta wing to the 


( 3 ) 


The elliptic integral E*(BC) depends only on the par ame ter BC ■ 

tan p 
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tangent of the Mach angle) and la therefore a constant for a given ving 


E"(BC) 


E* (BC) 


"a 


m* 

a 


DERIVATIVES 


■ / m_ 


s in reference 1 for C, t and C m . , 

"a a, 

0, thus merely acquire a factor 

■ E"(BC) 

(4) 

t _ 5| E n (BC) 

(5) 

, C T , and 
L \ 

\ 

and C, are 

* 2 

derived in reference 5. 

-ff O(B0) 

(6) 

~ H(BC) 

(7) 

- Si i(bc) 

32 

(8) 


where 


G(BG) 


1 - B^ 2 


(1 - SB^JE^BC) + B^SF^BC) 


H(BC ) = 3G(BC) - 2E" (BC ) 

2(1 - B^ 2 ) 


I(BC) = 


(2 - B^^E^BC) - BSc^f^BC) 


(9) 

( 10 ) 

( 11 ) 
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and F*(BC) and E * (BC ) are the complete elliptic intefgralaof the first 
and second kind s, respectively, with modulus k = ^1 - B2C2. 


DERIVATIVE C 

The pressure distribution over a thin delta wing in yaw (sideslip) at 
an angle of attack at supersonic speed has been obtained in reference 6 
and unpublished work. If the angle of yaw is assumed to be small ip « l? 

the rolling-moment coefficient can be expressed in the approximate form ^ 


C ~ E"(BC ) 

3 

Thus, the derivative with respect to p is 


c, = - E” (BC ) (12) 

*0 3 

An alternative derivation based on the surface potential, equation (3), 
for the unyawed wing will be given because the method provides the starting 

point for a derivation of C 7 , C_ , and C v . 

A r “r 

The potential for the disturbance velocity may be expressed relative 
to axes alined with the stream (wind axes) or with respect to axes that 
yaw with the body (body axes). For small angles of yaw (p « i\, the 

linearized equation for the potential has the same form relative to either 
system of sixes. The potential is determined by the normal velocity of 
points of the surface and by the orientation of the surface j for 
negligible thickness, this normal velocity is Just aV for all angles 
of yaw. The potential expressed relative to wind axes thus varies as 
the wing yews relative to these axes. The potential expressed relative 
to body axes is constant for small yaw because the orientation of the wing 
relative to the axes does not change. 

For wind axes, Bemoulli , s law has the form 

&P = 2pV M 
ox 

and the change in the pressure distribution with yaw results from the chan^p 
in the potential function with yaw. For body axes with small yaw, 
Bernoulli^ law has the approximate form 
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AP = 2pV(2£ - 0 

vdx ay/ 


(13) 


and tiie change In pressure distribution with yaw results 


term 



since 


0 does not change. 


from the 


In reference 1 in the section entitled "DERIVATIVE C 7 the 

0 

derivation employs body axes and equation ( 13 ) of the present paper. 

The surface potential used (equation ( 3 ) of reference 1) is the 
approximation for narrow vertex angle. Equation ( 3 ) herein for a 
general vertex angle may be used instead. Equation ( 3 ) herein differs 
only in the factor l/E*(BC), whence the earlier expression for C 7 

(equation ( 19 ), reference 1, with r = 0°) acquires this factor to agree 
with equation (12). 


DERIVATIVE C 

1 r 

The foregoing discussion of the triangular wing in yaw (sideslip) 
may be extended to provide a preliminary treatment of the case of a 
small angular velocity of yaw r. The corresponding extension for 
narrow vertex angle is made in reference 1. The treatment is general- 
ized to an arbitrary vertex angle for supersonic speeds, as before, by 
using equation ( 3 ) herein for the surface potential. Two changes then 
appear in the pressure equation, equation (20), of reference 1. The 
right-hand side is divided by E* (BC), and the term aC = xC 2 must be 
retained, since C 2 is no longer small compared with unity (C = tangent 
of semi vertex angle). With these changes, the derivation leads to 

C, = jta (— + ~Ae"(BC) (14) 

l r \9A 16 J 

In the derivation of equation (14), the spanwise variation in local 
Mach number caused by yawing is not taken into account although the 
variation in forward speed is taken into account. The surface potential 
that is used, equation ( 3 ), satisfies the linearized equation for a flow 
of uniform Mach number. This potential is inadequate to describe the 
compressibility effects associated with a spanwise variation of Mach 
number. 

Thus, consider a high— aspect— ratio rectangular wing with tips cut 
off along the Mach lines. In straight flight the Ackeret theory can be 
applied. The pressure difference is given by 
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AP = a 


2 pM 2 (Speed of sound ) 2 
\l M 2 - 1 


(15) 


In yawing flight the forward Telocity varies linearly along the span. If 
the rate of yaw Is made sufficiently low, the variation from wing tip to 
wing tip can he made so small that the flow is still nearly two-dimensional 
at any point. Thus the Ackeret theory is still applicable if the local 
Mach number is used at each s panwise station. 

The variation in pressure with local Mach number can be obtained from 
equation (15). As the Mach number is increased, the pressure decreases 
from infinity at M « 1 to a minimum at M = 1.4 and then increases again. 
Thus below Mach number 1.4 the faster moving sections of the yawing wing 
have the lesser lift. This result is contrary to subsonic behavior and to 
that which would be predicted if the spanwlse variation of Mach number were 
neglected. Thus the spanwise variation of the compressibility effect causes 
a reversal of the sign of the rolling moment due to yawing for rectangular 
wings at Mach numbers between 1 and 1.4, and at M = 1.4 the moment is 
zero. (This re stilt refers to yawing In a system of stability axes, fig. 3. 
For body axes, fig. 2, the effect is similar but the reversal extends 
to M = °°. ) 

A yawing triangular wing may be expected likewise to show an effect 
of the spanwlse variation in Mach number. If the triangle Is contained 
within the Mach cone from its vertex (the only case considered in this 
paper), however, the effect should be very much less than for the rec- 
tangular wing. In particular, where the predicted effect for the rec- 
tangular wing is a reversal of the sign of the rolling moment, the effect 
for the triangular wing Is expected to be merely a change in the magnitude. 

A reversal in sign is not expected until the edges of the triangle protrude 
from the Mach cane. This behavior Is inferred from the fact that the 
analyses of references 2 to 7 show many subsonic characteristics for 
triangles within the Mach cone and a marked change in characteristics for 
triangles with side edges outside the Mach cone. 


DERIVATIVES C v and C„ 
r p 


Extensive changes are necessary to generalize the treatment of C Y 

P 

and S In reference 1 to arbitrary vertex angles for supersonic speeds; 
therefore, the revised derivation is given in detail. 

The derivatives C y and C„ relative to body axes for a very 

P P 

thin delta wing without dihedral arise entirely from suction on the wing 
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side edges. Consider a condition for which the induced velocity normal 
to the edge is of the form 


▼» = 



(16) 


in the immediate neighborhood of the edge, vhere s is the perpendicular 
distance from the edge and K is a constant. Beference 3 points out 
that for such a flow there is a suction force per unit length of edge. 


f = 


it 


pK 2 \[l - M » 2 


(17) 


so long as the delta wing does not protrude from the Mach cone from its 
vertex. In equation (17), M* is the Mach number of the c omponent of 
the stream flow normal to the leading edge. The radical \J 1 — M*2 is 
the Prandtl-Glauert compressibility factor for the normal component of 
flow. Equation ( 17 ) is limited to real values of the radical by the 
condition expressed for the Mach cone. 


For the delta wing in rolling motion the induced velocity component 
has been obtained in reference 5 as 


U 1 



Angle of attack gives the additional contribution (reference 2) 

oVC 2 

*2 = + 


XU = ± 


E*(BC) ^ C 2 - 


The total iniuced velocity on the upper surface is thus the sum of 
anil Ug with the plus sign 



Very near the side edge this velocity is approximately 
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u 



—22— ± £2* i(bc) 
E*(BC) 2 


where the plus sign refers to the right edge and the minus sign to the left 
edge. 

If a similar calculation is made for t = it is found that as 

the side edge is approached the resultant induced Telocity yu 2 + t 2 
becomes normal to the edge. Thus the normal Telocity near the edge is 


\Ji7 


u 


The perpendicular distance of point (x,y) from the side edge is 

..ffcdi 

\Ji + cS 


The resultant induced Telocity Tery near the edge may therefore be 
expressed approximately as 



aV + l(BC)pCx 
E*(BC) ” 2 



which is of the form of equation (l6). The suction force per unit length 
of edge is from equation ( 17 ) thus 


f = 


pCx < 


a2yg + I(BC)_ p 2 C 2 x 2 + l(BC)aYpCx 
[e^bc)] 2 ^ ” E*(BC) 


[ \J( 1 + C 2 )(l - M* 2 ) 



the minus sign refers to 
can be reduced to 


The lateral component of this suction force is giTen by 
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1 - fj ( f E - f l) «* 

o a I(BC) vA - B^ 2 
= | pC 2 c 3 oVp — i-V 


E*(BC) 


The lateral-force coefficient is formed "by division "by - p^S, and 

2 

the derivative with respect to p'b/2V‘ is the stability derivative Cy 
It is p 


_ 2 m X(B0) \Jl - B 2 C £ 


e*(bc) 


(18) 


The yawing moment of the leading-edge suction about the vertex of 
the triangle is 


+ C 2 ‘dfx \ll + C 2 


- - S ♦ o 2 ) 


The moment about the reference point 


N = N-. + — cY 




is 


. - «. ^Wl + 9C 2 ) 

36 E* (BC) 


The yawing-moment coefficient is formed by division by ^ pV^Sb, and 

the derivative with respect to pb/2V is the stability derivative C- 
It is 




0n _ :-Jh- * *3) i mJIIj iE: 

HP \9A 16/ e*(BC) 


( 19 ) 



Ik 
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DERIVATIVES 


/ 


According to the discussion' on 


8 Ides lip' does not/alter the surf ad 

relative to hody axes. As d consequence the induced velocity distrlbu— ( 
t ion' is unchanged. Thus jrfie initially symmetric distribution of leading— 

P / 
rdd" i 



and C 

smal l angle of yaw or 
potential expressed 


\ 




and C 


n 


edge suction persists ip' sideslip, and the derivatives 
a:j*e zero. / \ 

( The surface potential relative to body axes is l^fcewise unaltered in 
the first approximation by a small, angular velocity of yaw. ^Accordingly, 
insofar as the pressure forces are concerned, the derivatives Cy and C 

are zero. Subsonic experience, however, suggests an appreciable 
Cj^-derivativer (damping In yaw) from profile drag. This damping derivative 

has been evaluated in reference 1 as 


“r 


\ 


\ 


C„ * 




\ 


V 

\ 


( 20 ) 


RESULTS AND DISCUSSION 


The formulas that have been obtained for the various stability 
derivatives are collected in table I. Derivatives obtained elsewhere 
are included for completeness, and the source is indicated in each 
instance. Expressions are given for two systems of coordinate axes. 

In the first column are shown the derivatives relative to the principal 

p 

body axes of figure 2 with origin a distance from the vertex of the 

triangle. In the second column are shown the results relative to 
stability axes with origin a distance x cg ahead of the ^c point. 

The relationship between the two systems of axes Is shown in figure 3. 
Equations for transforming from body axes to stability axes are given 
in reference 8; the shift in origin results in additional terms. 

In the transformation of the present results from principal body 
axe 8 to stability axes terms of order A ^/16 * and the more important 
terms of order a P are retained (see footnote, table I), whereas in 
reference 1 such terms are dropped as a consequence of the narrow vertex- 
angle approximation. - 
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These results for an arbitrary vertex angle may be compared with the 
asymptotic values for the case of vertex angle approaching zero given in 
reference 1. The present results for principal axes are found to differ 
from the asymptotic values (except for small terms in A^ and a^) only 
In the acquisition of certain factors which are functions of BC. Thus 
the asymptotic values for Ct , Ct . , C,, j C, , and C 7 are multiplied 

■ Lj a q 4 3 4 r 

by E"(BC) ; CL is multiplied by G(BC); Ct is multiplied by H(BC); 

q q 

C, is multiplied by l(BC); and C_ and C y are multiplied 
4 P P P 

by ■ = j(BC) . The parameter BC = s is the ratio of 

E* (BC) tan » 

tangent of the semivertex angle of the triangle to the tangent of the 
Mach angle. BC approaches zero, therefore, as the vertex an gle approaches 
zero. The several functions E"(BC), . . ., J(BC) all approach unity 
as BC approaches zero, and thus the derivatives obtained herein approach 
the asymptotic values of reference 1 as the vertex angle goes to zero. 

The variation of these stability derivatives with Mach number is 
contained entirely in the factors E” (BC), . . ., J(BC). The five 

factors are plotted against BC = ^ ■ an — » , the ratio of the tangent of 

tan u 

the semi vertex angle to the tangent of the Mach angle, in figure 4. 

The derivatives apply to a wing of triangular plan form and zero 
thicimess. The calculations are based on the assumption of potential 
flow with small disturbances, except in the case of the derivative C n , 

in which skin friction is considered. The predicted infinite negative 
pressure acting on an edge of zero thickness to yield a finite suction 
force is, of course, a mathematical idealization. (The local violation 
of the assumption of small disturbances is not serious.) Subsonic 
experience indicates that with a suitably rounded edge a considerable 
leading— edge suction force may be realized in practice, with the 
theoretical value an upper limit. On the other hand, a sharp leading edge 
is known to cause loss of the leading— edge suction. The requirements of 
extreme thinness and a rounded leading edge (that is, appreciable radius 
of curvature) are evidently in conflict. Thus the degree of applicability 
of the yawing-moment and lateral— force derivatives to actual triangular 
wings is 'uncertain. A further limitation on validity, already elaborated 

on in the section on C^ , exists also for the derivatives with respect 

t r 

to yawing velocity. The analysis neglects the spanwise variation in 
Mach number caused by the yawing (but not the spanwise variation in 
velocity). The result is an error in the magnitude of the yawing 
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derivatives that is expected to vary from zero for BC—^O to an important 
amount for BC — 3>1. 

The potential 0 satisfies the linearized equation of motion for 
the steady state hut not the more general linearized equation for unsteady 
motion. This circumstance implies that the present expressions for the 
stability derivatives are suitable only for steady motions, motions with 
small accelerations, or sinuous motions of low frequency. This limitation 
is accepted in all stability work and may become serious only in cases of 
high— frequency oscillations such as flutter. 


Langley Memorial Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va. , November 6, 19^7 
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1 ln the transformation from body axes t«nu of order a 2 hare teen neglected In comparison with unity, but tens* of 
order a 2 / A here teen retained, since they may be appreciable for em ail values of A. — •pp*' 
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Figure 1.- Axes and notation used in analysis. 



Figure 2.- Velocities, forces, and moments relative to 

2 

principal axes with origin at jc . 



Figure 5 Velocities, forces, and moments relative to 

2 

stability axes with origin at — c - x c ~. Principal 

5 s 

axes of figure 2 dotted in for comparison. 
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Figure 4. _ Elliptic tn fegra / factors of the 
stability denvaf/ves that determ me 
their var/af/on with Mach number. 

( See Fable X. ) 





